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ABSTRACT 

Transporting solids of different sizes is an essential process in the evolution of protoplanetary 
disks and planet formation. Large solids are supposed to drift inward; high-temperature minerals 
found in comets are assumed to have been transported outward. From low-gravity experiments 
on parabolic flights we studied the light-induced erosion of dusty bodies caused by a solid-state 
greenhouse effect and photophoresis within a dust bed's upper layers. Tje gravity levels studied 
were 0.16c/, 0.38g, lg, and 1.7 g. The light flux during the experiments was 12 ± 2 kW m~ 2 and 
the ambient pressure was 6 ± 0.9 mbar. Light-induced erosion is strongly gravity dependent, 
which is in agreement with a developed model. In particular for small dusty bodies ((sub)- 
planetesimals), efficient erosion is possible at the optically thin inner edges of protoplanetary 
disks. Light-induced erosion prevents significant parts of a larger body from moving too close 
to the host star and be being subsequently accreted. The small dust produced continues to be 
subject to photophoresis and is partially transported upward and outward over the surface of the 
disk; the resulting small dust particles observed over the disk's lifetime. The fraction of eroded 
dust participates in subsequent cycles of growth during planetesimal formation. Another fraction 
of dust might be collected by a body of planetary size if this body is already present close to the 
disk edge. Either way, light induced erosion is an efficient recycling process in protoplanetary 
disks. 

Subject headings: methods: laboratory-, planets and satellites: formation - protoplanetary disks 



1. Introduction 

The formation of planets in protoplanetary 
disks is strongly tied to the existence of dust 
which m akes up abo ut 1% of the total mass of 



the disk (Nat ta et al.1 120071 ). The basic collisional 



growth model assumes that, as a first step toward 
planet formation, kilometer-size planetesimals are 
formed by sticking collisions of dust particles 
which are only held t o gether by surface forces 
(iBlum fc Wurml 120081: IWeidenschilling fc Cuzzl 
19931 : IWindmark et al.1 l2012t ). There are some 
issues in planetesimal formation which are not 
yet understood. Somewhat problematic is inward 
drift especially of meter-size bodies. In a standard 
disk the pressure is supposed to decrease with dis- 
tance from the star ( Havashi et al. 1985). The gas 



is therefore supported by a pressure gradient and 
rotates with sub-Keplerian velocity. Small solids 
couple to the gas and thus also rotate more slowly 
than Keplerian. However, solid bodies are not 
supported by the pressure gradient. This leads to 
an inward drift of the solid bodies which - for a 
meter size body in a minimum m ass nebula - can 
be as large as 1 AU in 100 years (jWeidenschilling 
Il977l ). Such an inward mass transport of solids 
leads to a redistribution of matter: the outer part 
of the disk becomes depleted in solids, while the 
inner part becomes enriched. The physics of the 
inner 1 AU of protoplanetary disks strongly de- 
pends on the fate of this incoming matter: it 
might contribute to the local planet formation 
process and add mass to forming planetary bod- 
ies, it might evaporate and be accreted by the star, 
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or it might be recycled in some more complex way. 
As the physics in the inner part of the disk is com- 
plex the evolution of solid bodies is still an open 
question. To address this problem, we started 
to investigate a mechanism of particle recycling 
which is especially active in this inner region of a 
protoplanctary disk and, more specifically, close 
to the inner edge of a disk within the inner 1 
AU. There is a growing number of observations 
which show that, in some phases, protoplanetary 
disks contain inner gaps which are optically thin 
but still contain gas, while the outer part i s still 



dense and optical l y thick ( D'Alessi o et al 



C alve t et al.l 120021: ICurrie fc Sicilia-Aguilar 



2005 



2011 



Sicilia-Aguilar et all 20081) . In other words, solids 



at the inner edge of a disk are embedded in a 
gaseous environment and are illuminated by the 
star. 



Wurm fe Kraussl (|2006l) first described a mech- 
anism that significantly erodes dusty bodies to 
small particles under very general conditions, e.g. 
by mere illumination and at low gas pressu re (kW 



of intensity and mbar gas pressure). IWu 



(|2007l ) showed that this erosion process also works 
under the conditions of transitional disks. There 
will most likely be a small shell at the inner edge 
of a disk where the co nditions for particle erosion 
are met ( Wurml 2007 ). Every dusty body which 
enters this zone and which is too small to hold 
its dust by gravity (sub-planetesimal) is subject 
to erosion by this process. Different scenarios are 
possible as to how this contributes to the evolu- 
tion of the disk. The dust might, for example, 
be transported upward and outward over the disk 
by photo phoresis. Such a mode l has been sug- 
gested by IWurm fc Haackl (|2009h . At later times 
or throughout the optical thin inner parts the 
dust might also be transported by photophore- 
sis in the midplane (Wurml 20071: iTakeuchi fc Lin 
2003[ iHerrmann fc Kriyovl 120071: iMoudens eTal 
201ll:lMousis et al.l2007l ). iKelling fc Wurml (|20lTb 
showed that the erosion process is a suitable mech- 
anism to explain the existence of small dust par- 
ticles which are observed in protoplanetary disks 
over their entire lifetime. The dust might also 
be recycled locally and added - probably more 
efficiently than by the original larger bodies - 
to existing large planetesimals or protoplanetary 
bodies to boost their growth. These aspects are 



speculative and not the main focus of this paper 
which reports on experiments that analyze the 
gravity dependence of the light induced erosion 
mechanism providing basic input for future mod- 
eling. 

A sketch of the light induced erosion mech- 
anism is shown in Figure [TJ the illumination 
of a dusty surface first leads to a solid-state 

greenhouse (SSG) effect. This effect is known 
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to exist for ice and comets (Niedcrdorfer 193 



Davidsson fc Skorov 2002: Kaufmann et al. 200 



It means that visible light enters a porous body 
and is absorbed over some depth, but as the body 
is opaque to infrared radiation only the top sur- 
face can cool by thermal emission. Therefore, the 
maximum temperature is below the surface. In 
recent years it was sho wn that this is also t r ue for 



purely dusty surfac es (Kocifai et al.l 2010L 12011 



Kelling et all 120111 IKelling fc Wurml l201lh . If 



temperature gradients across a particle exist as in 
the case of an SSG profile, and if the particle is em- 
bedded in a gaseous environment, a photophoretic 
force acts on the particle from the warm to the 
cold side. Hence, the surface particles of an illumi- 
nated dust bed are subject to a lifting force. If this 
force is stronger than gravity and cohesion, aggre- 
gates are ejected. The SSG and photophoresis are 
the basis for light induced erosion of dusty bod- 
ies. The physics be hind this has been studied in 
a number of works ( Wurm fc Krauss 20061: Wurm 



[ 20071: IKelling and Wurmll2009l: IKelling et alj|20ll 
IKelling fc Wurmll201ll ). 



In protoplanetary disks, dust particles on dusty 
sub-planetesimals are bound by cohesion and grav- 
ity. However, the gravitational force is small due 
to the small mass of the objects. Therefore pho- 
tophoretic forces acting on the particles have to 
overcome less gravity and the effect of dust erosion 
induced by illumination is more intense compared 
to larger bodies or laboratory experiments. In the 
first microgravity experiments Wurm et al. (2008) 
showed that the threshold light flux for parti- 
cle erosion indeed depends on the gravity level. 
At low gravity cohesion eventually becomes the 
dominant opposing force. The erosion rate was 



not m easured in the experiments of IWurm et al 



(|2008l ): however, this is important in determining 
how rapidly a body can be eroded or what mass 
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flux it might provide in small particles if it crosses 
the erosion zone in a protoplanetary disk. 

To quantify the effect of light induced erosion 
for application to (sub)-planetesimal recycling we 
carried out low gravity experiments and deter- 
mined the gravity dependence of the erosion rate. 
The experiments were carried out during the first 
Joint European Partial- G Parabolic Flight cam- 
paign. On this campaign, over 3 days, 13 parabo- 
las at 0.1 6g, 12 parabolas at 0.38g and 6 parabolas 
at Og (g — 9.81 m s -2 ) were provided. We also 
developed a preliminary model which agrees well 
with the experimental results. 

2. Experiment 

A sketch of the experimental setup is shown 
in Figure [21 Four separate chambers are placed 
in one vacuum chamber. Each chamber contains 
a dust sample with a depth of 2 cm covering a 
total area of 6 cm 2 . For the campaign reported 
here the chambers were filled with basalt pow- 
der with a broad size distribution between and 
125 /im. We consider this dust particle size range 
suitable to simulate the dust in protoplanetary 
disks. They incl ude micron-size particles discussed 
Brearlevl (|l999t ) with reference to the matrix in 
chondrites (typically a few microns but q ccasion- 
ally up to about 20 fim (jBrearlev! fl999h ). The 
constituents in interplanetary dust particles or 
com etary material are smaller, on the order of 0.1 
/um (jWozniakiewicz et al.l 120121 ) which might not 
be adequatley represented in our experiment sam- 
ple. On the other hand, chondrites largely consist 
of chondrules of submillimeter size, so the larger 
grain fraction in our sample might correlate with 
these dimensions. Certainly, the size dependence 
of the erosion mechanism should be considered in 
more detail in the future. 

The dust samples were dried at about 500 K for 
24 hr and were stored in a desiccator until they 
were placed in the vacuum chamber just before 
the experiments. The ambient pressure of air in 
the vacuum chamber was then set to p = 6 ± 0.9 
mbar for all experiments. 

During each experimental run one sample was 
illuminated by a diode laser of 2 W optical output 
at 655 ± 10 nm. The laser profile was homoge- 



nized by coupling the laser to an optical fiber and 
projecting the fiber outlet onto the sample surface 
with a size of 25 ± 4 mm 2 . The average inten- 
sity was 12 kW m~ 2 which varied by 17 % within 
the spot. This corresponds to the inner regions of 
a protoplanetary disk < 0.4 AU where the light 
flux for a star similar to the sun is I > 10 kW 
m~ 2 . An accelerometer was used to measure the 
residual gravity acting on the experiment. The ex- 
perimental chamber was supported by vibrational 
attenuators to prevent any influence of airplane 
vibrations, which might reduce the cohesion be- 
tween the dust particles. Released aggregates were 
observed from the side by a camera at 60 images 
per second. 

Measurements (images) were taken for lg, the 
acceleration phase at 1.7 g and the low-gravity 
phase at Q.lQg (Moon level) or 0.38g (Mars level). 
To confine the dust within the individual cham- 
bers, they were enclosed in a glass cylinder. This 
cylinder and the chambers were independently 
rotatable by a rotational feedthrough from the 
outside. After each parabola, and hence before 
every new experiment, the inside of the rotating 
glass cylinder was cleaned by a small strip of foam. 
This ensured that the laser and the camera view 
were not affected by dust. Rotating the inner, 
dust containig part of the chambers ensures that 
for every new experiment a fresh spot of the dust 
surface was illuminated. This prevented selection 
effects, e.g., that all ejectable particles (e.g. of a 
certain size fraction) were already released. Af- 
ter removing the background noise of the camera 
images the average gray scale value was taken as 
a relative measurement for the ejected mass or 
particle number. This assumed that the ejected 
particle distribution stays constant with time, that 
averaging is done over a significantly large num- 
ber of particles and that the image is not satu- 
rated with particles. As we provide new spots 
for each experiment, observe well separated ag- 
gregates and take an average over a large number 
of aggregates in typically 900 images for a given 
g-phase (Figure [3]), we consider the gray value as 
a suitable measure. Figure [3] shows an overlay of 
100 images of the ejected aggregates at different 
g-levels. The number of ejected aggregates clearly 
increases with lower gravitational acceleration. 
A parabolic flight has the advantage that dif- 
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Fig. 2. — Experimental setup: four chambers filled with dust samples are enclosed by a rotatable glass 
cylinder in a vacuum chamber. A homogeneous laser spot of 25 ± 4 mm 2 size illuminates the dust bed with 
an intensity of 12 ± 2 kW m -2 . The particles being ejected at different ^-levels are observed by a side view 
camera. 
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Fig. 1. — Visible radiation (laser) enters the dust 
bed and leads to a solid state greenhouse (SSG) 
effect, with a maximum temperature below the 
surface (indicated by the black line: larger tem- 
peratures are represented by thicker lines). Pho- 
tophoresis F p h acts on the surface particles. When 
this force overcomes gravity and cohesion, the sur- 
face particles are lifted. A crater is formed where 
again temperature gradients develop. If the illumi- 
nated spot is limited in size, the developing crater 
might prevent further release of particle. 
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Fig. 3. — Dust erosion of basalt powder at 
6 mbar ambient pressure and a light intensity of 
12 kW m -2 at four different gravitational accel- 
erations. Each figure is an overlay of 100 images 
to emphasize the increase of particle release with 
lower accelerations. 
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ferent gravity levels can be tested with the same 
dust sample in a short period of time. However, 
the effects of the transitions between the differ- 
ent g-levels have to be filtered out. They induce 
a higher erosion rate during the transitions from 
high-g to low-g. This is due to the release of grav- 
itational and photophoretic tension from the dust 
bed, as higher gravity levels require a steeper tem- 
perature gradient for erosion. It takes a certain 
time before the new equilibrium temperature pro- 
file is established at which aggregates are ejected 
regularly. Figure |4]shows which data were selected 
to determine the number of ejected aggregates in 
the different phases. In the lg-phase of every 
parabola the laser and camera were turned on. 
First erosions occur after 2-3 s. During the tran- 
sition from 1.7g to lower ^-levels, massive particle 
releases appear (peaks in Figured]). These releases 
are an artifact of the tension release as previously 
mentioned and are not used to determine the ero- 
sion rate. After the transition particle releases 
occur continuously at lower g-levels for about 10 
s. The dust bed surface changes as the laser spot 
removes particles and creates an erosion crater 
(see also Figure [lj . The presence of a crater rim 
induces edge effects, which efficiently reduce the 
erosion rate, especially in the low gravity phase 
with high erosion rates. Therefore, we do not in- 
clude the later times for averaging. After the data 
reduction, a measure of the erosion rate is given 
as the average brightness of all images within a 
single g-phase. As the brightness is an arbitrary 
measure of the erosion rate, the different ^-phases 
of one parabola are scaled to a value of 1 at lg. 
The parabolic flight provided results for 12 
parabolas: 6 Moon parabolas and 6 Mars parabo- 
las. Figures [5] and 6 [5] show the data sets for the 
Moon and Mars parabolas. 

We further average over all parabolas (6 data 
points for 0.16<? and 0.38<? and 12 data points 
for lg and 1.7g) and obtain the average ejected 
mass rate over gravity. In Figure [SJ all (averaged) 
experimental data are presented in the context of 
the model developed below. 

3. The Model 

The continuous aggregate ejections from the 
dust bed are caused by an inverse temperature 
gradient below the dust bed's surface and a 
photophoretic force F p h (|Wurm fc Krauss 2006; 




Fig. 4. — Example of the collected raw data of a 
Moon parabola. The black line shows the grav- 
itational acceleration in units of g = 9.81 m s~ 2 . 
The gray peaks indicate the eroded mass in ar- 
bitrary units. For the analysis only the data in 
the green boxes were taken where the release of 
aggregates is continuous. 
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Fig. 5. — Dataset for the 6 Moon parabolas. The 
particle erosion rate is set to 1 at lg. The results 
from the experiments are presented in gray, the 
average value in black and the black line shows 
the trend for the averaged data according to the 
model developed in Section 3. 
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Kelling et al.ll201 libelling fc WurnfeOllI) . In de- 
tail the idea behind the lig ht induced ag g regate 



Kocifai et all (j2010L 



eject ions is the following: 
20111 ) showed, that if a dust bed is illuminated 
from the top, the maximum temperature T max is 
up to some 100 ^m below the surface - this is 
called SSG. Hence, there is a temperature gradi- 
ent dT/dz pointing from T max toward the cooler 
surface (inverse temperature gradient; z as depth 
within the dust bed). Typical thermal conduc- 
tivities Kd of the dust bed are on the or der of 
n d ~ 10~ 2 W m^ 1 K" 1 (|Krause et al.ll201ll) . The 
inverse temperature gradient then takes values of 
dT/dz ~ 10 5 K m" 1 at 10 kW m~ 2 illumination 



([Kocifai et al.ll2011l) . 

The upper most aggregates have a temperature 
gradient along their surface. Small particles with 
a temperature gradient in a gaseous environment 
are affected by a photophoretic forc e which can 
be w ritten for spherical particles as (jRohatschek 
19951) 



2F„ 



+ 



F, 



= Dj-a 2 — 

oz 



2 

dT 



(1) 

(2) 
(3) 



where p is gas pressure. F max and p max are par- 
ticle and gas dependent parameters determining 
the pressure p max at which the photophoretic force 
-Fmax is a maximum. D is a gas dependent param- 
eter, a = 1 is the accommodation coefficient, a is 
the particle radius and dT/dz is the temperature 
gradient over the particle's surface. Particles in 
the /im regime are most strongly affected by pho- 
tophoresis at mbar pressure. 

Photophoresis acts on the upper-most aggre- 
gates in illuminated dust beds. The threshold 
for particle ejectio n can be written as follows 
(jWurm et al.ll2008l ): 



^ P h th = Am ■ g c 



(4) 



with gravity Fq = Am • g cn v, where g env is the en- 
vironmental gravitational acceleration, Am is the 
ejected mass, and Fc is the cohesion force. If the 
photophoretic force F p h overcomes gravity and co- 
hesion, aggregates are released from the dust bed's 
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Fig. 6. — Dataset for the 6 Mars parabolas. The 
particle erosion rate is set to 1 at lg. The results 
from the experiments are presented in gray, the 
average value in black and the black line shows 
the trend for the averaged data according to the 
model developed in Section 3. 
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Fig. 7. — Trajectories of dust aggregates leaving 
the surface of an illuminated dust bed (basalt, 0- 
125 fj,m). The experiment was carried out on the 
ground at 6 mbar ambient pressure and at an in- 
tensity of 12 kW m -2 . 
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surface. The mass ejection rate N = Am/ At with 
At as time interval depends on the time scale at 
which the threshold temperature gradient (see e.g. 
Equation fl2J)) within the illuminated dust bed is 
reset after an ejection. The temperature evolution 
toward equilibrium withi n the dust bed follow s an 
exponential dependence (jKocifaj et al.l 120111 1. If 



the threshold gradient to induce particle release is 
reached after a time span At much smaller than 
the exponential timescale (large light flux) the 
temperature gradient can be approximated by a 
linear increase with time 



dr. . 

oz 



L0 2 t. 



(5) 



Therefore the threshold temperature gradient af- 
ter a time span At can be written as 



— =— {At )=uj 2 At. 

OZ th OZ 



(6) 



With Equation ([5J the photophoretic force (Equ- 
tation j3j) can be written as 



^ P h th = w 3 Ai 



(7) 



with UJ3 = uj\ijj2- Hence, after the time interval At 
the threshold photophoretic force (F p h th ), Eq.dU) 
is reached 



-Fph th = uj 3 At — Am ■ g cnv + Fc, 



(8) 



and aggregates are released in steps always after 
an average time 



At 



Am ■ g CDV + F c 
w 3 



The mass ejection rate therefore is 

w 3 



AT/ . Am 

N(genv) = 



<7env + ^„ 



(9) 



(10) 



Equation (|10[) was used to fit the experimental 
data, and gives W3 = 1.6 m s~ 3 and Fc / Am = 6.1 
± 1.2 m s~ 2 . u>3 is scalable arbitrarily as it ac- 
counts for the nature of the brightness measure- 
ment as the relative measurement for one part. 
However, Fc / Am contains information on the co- 
hesive force by which an ejected mass fraction is 
bound to the dust bed. 

Due to the low resolution of the analyzed images 
we can not give an accurate assumption for the 



eroded particle mass Am. Hence, we can not esti- 
mate the exact value of Fc from our experiments. 
The acceleration Fc/ Am from the model fits the 
accelerations measured in first follow-up labora- 
tory experiments, shown in Figure [7] basalt pow- 
der with a broad size distribution between and 
125 /im was placed in a vacuum chamber with 6 
±0.9 mbar ambient pressure. It was then illumi- 
nated by a laser with 12 ± 1 kWm~ 2 and a spot 
size of 5 ± 0.5 mm radius. The resulting ejection 
of dust aggregates was measured at a camera with 
500 images per second, so the accelerated aggre- 
gates could be tracked. We measured accelerations 
of 14.4 ± 4.9 m s~ 2 , averaged over nine ejected 
aggregates, which is only twice as high as the ac- 
celeration measured on the parabolic flight. This 
factor two can be explained by the fact that we 
could not use exactly the same dust sample that 
we had on board for the parabolic flight. 
To visualize the influence of the cohesion we plot 
the ejection rate of Equation (fTU|) divided by the 
factor 0J3 / (Fc / Am) for the experimental data (av- 
erage values of the ^-phases) which effectively nor- 
malizes the experimental data to 1 at Og. In Figure 
[HI the value Fc / Am is then varied by a factor / to 
show the general dependence on Fc / Am. There is 
a prominent increase in the ejection rate towards 
smaller Fc / Am where it becomes more and more 
dominant at g cm < I5 while cohesion is negligible 
for g cnv > lg. 

This simple model neglects self-shielding of the 
dust bed by ejected aggregates. The more dust 
that is ejected, the more light is shielded from 
reaching the surface. Depending on the relation 
between the ejection rate and the rate at which 
aggregates are transported away from the ejection 
location, e.g. by gas drag in protoplanetary disks, 
the optical thickness of the layer above the dust 
bed's surface will increase. From the experiments 
and based on the assumption that the light flux 
might be reduced to less than 1/e of it's inten- 
sity, ejection rates that are increased by a factor 
10 likely still do not affect the ejection rate. This 
might set an upper limit to the ejection rate but 
it is only a rough estimate here which has to be 
quantified in the future. For planetary bodies the 
effect might lead to a size sorting over time where 
large particles fall back down while small grains 
get entrained into the atmosphere. 



7 



4. Applications and Conclusion 



The effect of light induced dust ejections can be 
a key process of recycling dust during the forma- 
tion of planets in protoplanetary disks. Planetesi- 
mals or smaller d usty bodies a t the inner edge can 
be partly eroded (|Wurml 120071 ). 

The position of the edge in a number of tran- 
sition, pre-t ransition and f ull di sks has been 
modeled by lEspaillat et al. ( 2012 ). Here, pre- 
transition disks relate to a disk with substantial 
matter within the gap modeled by an optically 
thin inner wall at the dust sublimation radius and 
a second outer wall. In Table [TJ we extract their 
modeled radii of the disks and the luminosities of 
the host stars. Based on the photophoresis the lift 
force on particles depends on the environmental 
parameters light flux /, gas pressure p and gas 
temperature T. In the free molecular flow regime 
(low pressure) Equation [1] simplifies and connects 
these three parameters to 



F 



(11) 



(|Wurml 120071 ). with a parameter b dependent on 
the dust bed. 

Our experiments show that particle eruptions con- 
tinuously occur at all gravity levels at the used 
light flux I = 12 kW m~ 2 , pressure p — 6 mbar 
and temperature T of about 300 K. We did not 
probe the threshold condition where no dust is 
lifted but our values are consistent with th reshold 
condit ions measured for eruptions at lg bv IWurm 
(|2007l ). He found I = 50 kW m~ 2 at 1 mbar and 
300 K. Therefore, we expect that light induced 
erosi on occurs as long as Ip/T > 15 kW m -2 Pa 
K" 1 (IWurmll2007l) . 



In a minimum mass solar nebula (H avashi et al 



19851 ) with given pressure, temperature and lu- 
minosity dependence this condition is reached at 
a distance o f about 0.4 AU or Mercury's orbit 
ll\VuriiJl2()07h. 

For observed extrasolar systems we only know 
the luminosities for certain, not the pressure or 
temperature at the midplanes. To determine the 
distance of the dividing line between erosion and 
stability for the different objects, we assume that 
the dividing line is where the light flux of the host 
star equals that of a solar-type star in a minimum 
mass solar nebula or 0.4 AU x y/L/L sun (Table[TJ. 
The implicit assumptions in this equation are that 
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Fig. 8. — Ejection rate normalized to 1 at Og for 
the experimental values using Equation 1101 The 
experimental value (squares) for Fc/Am is var- 
ied by a factor / with / = 2, 1,0.5,0.1,0.01 (bot- 
tom to top). The color gradient indicates that 
the more aggregates that are ejected, the more 
light is shielded from reaching the dust bed's sur- 
face. From the experiments, we estimate that an 
increase of the rate up to factor of 10 might be 
realistic for low cohesive powders. 
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Fig. 9. — An inward drifting dusty body (B) can 
be eroded in the optically thin part of a protoplan- 
etary disk by the effect of light induced ejection. 
The planetesimal loses mass which can be trans- 
ported outwards again by e.g. photophoresis. 
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Object Star's Luminosity (L sun ) Inner Wall (AU) Erosion Dividing Line (AU) Outer Wall (AU) 



FM 515 
FM 618 
LRLL 21 
LRLL 31 
LRLL 37 



2.5 
2.2 
3.8 
5.0 
1.3 



0.12 
0.22 
0.13 
0.32 
0.17 



0.63 
0.59 
0.80 
0.89 
0.46 



45 
11 
9 

14 

5 



Table 1: Pre-transitional Disks Modeled bv lEspaillat et al.l (|2012l ). The erosion limit distance is the distance 
up to which disassembly of larger dusty bodies would be possible. 



the gas pressure at this distance and the temper- 
ature are the same as in the minimum-mass solar 
nebula at 0.4 AU. 

For the temperature this is a consistent assump- 
tion as the blackbody radiation temperature 
which defines the temperature of the disk in the 
minimum-mass solar nebula would be the same at 
the same incoming light flux. The pressure is less 
certain but one might argue that more luminous 
stars have denser disks. We therefore assume the 
pressure to be equal to that in a minimum mass 
solar nebula at the dividing line. However, more 
complex models would not change the following 
general statement. 

We always find the erosion limit within the gap 
of the pre-transitional disks. In general, erosion is 
possible at the inner edges of full disks, although 
melting or sintering of particle ensembles might 
be an issue which has to be studied in further 
detail. For full disks and pre-transitional disks 
erosion might clearly provide the small dust ob- 
served over the lifetime of the disk. As the erosion 
limit distance is closer to the star than the outer 
wall in pre-transitional disks, particles released by 
the erosion mechanism can drift outward to this 
second wall by photoph oresis and enhance the 
particle density locall y (jKrauss fc Wurm 2005; 
Haack fc WuTrnl l2007t) . This will induce an en- 
hanced growth of dust aggregates. At the outer 
wall the growing bodies are safe from destruction 
unless they are scattered or drift inward again to 
the inner erosion zone during further evolution 
and recycling. 

The holes of transitional disks are all larger 
than the erosion limit distance. As they have a 
largely reduced density it is unclear to what extent 
erosion will work here. The erosion limit distance 
would s hift inward sig nificantly. Earlier calcula- 
tions by IWurm ( 2007 ) showed that, for example, 



in the case of TW Hya erosion close to the star 
would still be possible but this strongly depends 
on the gas content. The erosion zone likely de- 
creases in thickness with time. However, even if it 
were only a very thin zone, dusty bodies traversing 
it might be significantly eroded. In the optical but 
thin still slightly gaseous hole, the dust would be 
transported by photophoresis and locally increase 
the particle density again. 

Therefore, at early full disk times erosion will 
influence the particle recycling at the inner edge. 
In later pre-transitional and transitional evolution 
times of the clearing phase the light induced ero- 
sion will be an important factor in shaping the 
size distribution or even triggering later formation 
of larger bodies at the outward moving edge. 

The experiments reported here show that the 
number of ejected particles increases strongly with 
decreasing ^-levels. The experiments performed 
on a parabolic flight display an inverse dependence 
on gravity and specific cohesion. The ejection rate 
is higher by a factor two in 0.16<7 than in Ig for 
the basalt sample used. The importance of gravity 
deduced from these data can more clearly be seen 
in the model. If cohesion is reduced the ejection 
rate might be up to one or two orders of magnitude 
more efficient at low-<?. Lower cohesion does not 
change the ejection rate strongly at high g-levels. 

As a consequence of microgravity on small plan- 
etesimals or sub-planetesimals, the effect is more 
intense than observed from groun d-based experi- 
ments so far. iKelling et al.l (|201ll ) estimated from 
lg experiments a mass- loss rate of 10 -5 kg s _1 
m~ 2 or about 10 3 kg year -1 for an area 1 m 2 
in size. Low gravity conditions would boost this 
value by a factor of a few. This compares well to 
the mass of a porous dusty body of meter size with 
a density of 1 g cm" 3 . The total mass is some 10 3 
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kg. It is therefore possible to erode, in less than 
one year, a significant part of such a body if it 
consists of loosely bound dust close to the disk 
edge. As the erosion zone is a fraction of an AU 
(<Wurml[2007h the experiments show that an ob- 
ject that passes the edge will efficiently lose mass 
(Figure |9]). Even for objects that lose mass but 
are not completely destroyed the effect is impor- 
tant as it provides a non-gravitational force on the 
body and will change the orbital parameters; the 
details, however, are beyond the scope of this pa- 
per. 

The mechanism is more efficient in an atmo- 
sphere consistin g of hydroge n inste ad of air. The 
Equation [3] by iRohatschea (|l995l ) for the pho- 
tophoretic force shows a dependence on different 
parameters. The maximum force is proportional 
to a gas dependent parameter 
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This parameter contains the dynamic viscosity of 
the gas r) and the average thermal velocity c, which 
is proportional to the square root of l//x, the in- 
verse molar mass. Assuming ^ a ; r = 17.1 /iPa-s and 
rjH 2 = 8.4 /xPa-s, as well as ^ a j r = 29 g mol -1 and 
Atproto = 2.31 g/mol, we get D proto = 1.75 x D air . 
The maximum force is therefore 1.75 times higher 
in a protoplanetary disk than on Earth. 

In principle the erosion mechanism might com- 
pete with regrowth as a body moves through the 
disk and collects dust particles. This is an inter- 
esting aspect as it changes the surface morphol- 
ogy and the susceptibility for re-erosion. To what 
degree this process is important depends on the 
position of the dividing line with respect to the 
inner edge of the optically thick disk. If erosion 
only occurs in the transitional zone far away from 
the edge then this whole zone is optically thin by 
definition and re-growth is not important. This 
might be explained as follows. Assuming all dust 
particles to be spherical with a radius of 5 //m, 
the particle cross section is 78.5 x 10~ 12 m 2 . The 
disk would be optically thick if its cross section 
were completely covered with dust particles. At 
typical distances of 0.5 AU and with a scale height 
of 0.1 AU this would be 7 x 10 21 m 2 . Dividing 
by the particle cross section, this gives 9 x 10 31 
particles. Spread out over an estimated 0.5 AU in 
distance the particle density would be n = 0.17 



particles m~ 3 . The largest drift velocities con- 
sidered in disk models are on the order of 50 m 
s _1 . If a body of 1 m radius colleced dust at these 
speeds the growth rate would be 26.7 particle per 
second or with a particles density of 2000 kg m~ 3 
it would be 3 x 10 -11 kg s _1 . Compared to the 
erosion rate this worst case estimate shows that 
re-growth is negligible. 

This changes at the edge. If the transition zone 
is only, for example, 100,000 km thick the parti- 
cle density would be 127.3 particles m~ 3 and the 
mass growth 2 x 10~ 8 kg s" 1 . This is a factor of 
500 less than the erosion rate under perfect condi- 
tions, but for non-perfect conditions, e.g., further 
toward the dark side of the edge growth and ero- 
sion might be comparable. 

Another application of the effect at low gravity 
is found on Mars. By adjusting the conditions 
(light flux and temperature) to Mars, dust en- 
trainment in the a tmosphere might be explained 
(|Wurm et al. 20081) . The low-gravity experiments 
suggest higher erosion rates on Mars. 

Most dust samples that absorb visible radia- 
tion show light-induced ejection. We have quan- 
tified the gravity effect for a single sample so far. 
We have experimental evidence that cohesion in 
the context of this ejection mechanism can vary 
strongly for different dust samples, but quantifica- 
tion is subject to future research. Overall, erosion 
of a dusty body at the low-gravity conditions of 
small bodies at the inner edges of protoplanetary 
disks is a major recycling process of dust, and is 
important in the context of planet formation. 
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